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Abstract. We consider the evolution of a two-mode system of bosons under the 
action of a Hamiltonian that generates linear SU(2) transformations. The Hamil- 
tonian is generic in that it represents a host of entanglement mechanisms, which 
can thus be treated in a unified way. We start by solving the quantum dynamics 
analytically when the system is initially in a Fock state. We show how the two 
modes get entangled by evolution to produce a coherent superposition of vortex 
states in general, and a single vortex state under certain conditions. The degree 
of entanglement between the modes is measured by finding the explicit analytical 
dependence of the Von Neumann entropy on the system parameters. The reduced 
state of each mode is analyzed by means of its correlation function and spatial 
coherence function. Remarkably, our analysis is shown to be equally as valid for 
a variety of initial states that can be prepared from a two-mode Fock state via 
a unitary transformation and for which the results can be obtained by mere in- 
spection of the corresponding results for an initial Fock state. As an example, 
we consider a quantum vortex as the initial state and also find conditions for its 
revival and charge conjugation. While studying the evolution of the initial vortex 
state, we have encountered and explained an interesting situation in which the 
entropy of the system does not evolve whereas its wave function does. Although 
the modal concept has been used throughout the paper, it is important to note 
that the theory is equally applicable for a two-particle system in which each par- 
ticle is represented by its bosonic creation and annihilation operators. 

(Figures in this article are in colour only in the electronic version) 
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1. Introduction 

Nonclassical properties of quantum states are actively being studied for their relevance 
in quantum computation. It is known that quantum entanglement is the key to 
performing communication and information processing tasks that cannot be realized 
classically. For this reason, there has been a surge of activity towards preparing, 
identifying and quantifying entangled systems [Q. 

An important source of quantum entanglement has been the polarization- 
entangled two-photon states generated from type-II phase-matched parametric down 
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conversion 2 . A variety of other entangled states can be produced by using various 
polarizing components. More recently, the subject of quantum information processing 
has been given a new direction with the realization that a number of quantum logic 
operations can be performed using single photons and methods of linear optics 
Even a method for quantum teleportation was proposed and implemented 0] . Clearly 
one needs to examine, in full generality, the question of transformation of an arbitrary 
input state by a device which can mix different states. 

We note that a number of special cases for the generation of entanglement using 
linear optical devices have been investigated. Huang and Agarwal |S] considered 
multimode systems described by a Hamiltonian that is quadratic in the mode 
operators. They derived conditions for the generation of an entangled state when the 
input state was represented by a Gaussian density matrix. Their treatment covered a 
large class of states including squeezed coherent states and even states with thermal 
noise. However they did not consider the case of input fields in Fock states. More 
recently and more specifically, Kim et al [5] examined the question of the generation 
of entangled states by a beam splitter using Fock states as input fields. 

In this paper we specialise to intensity- or number-preserving linear transfor- 
mations belonging to the SU(2) group. For two-mode states characterized by the 
annihilation operators a and b, such transformations can be generated by evolution 
under a Hamiltonian of the form 

H = g(^be^ + h.c) +0(a t a-6 t 6) (1) 
where g and are real constants. Introducing the generators of the SU(2) group as 

J 1 = (a*b + atf)/2, J 2 = (a f 6 - atf)/2i, J 3 = {a) a - tfb)/2, (2) 
the Hamiltonian Q can be rewritten in the form 

H = vi Ji + V2J2 + V3J3 (3) 

with 

V\ = 2g cos 0, t>2 = — 2gsin</>, v% = 2Q. (4) 

Motivation for the present work comes from the realization that a Hamiltonian 
of the form Q can represent a host of entanglement mechanisms which can thus be 
treated in a unified way. Several examples are given as follows. 

The beam splitter, used by many authors as an entangler [7] can be described by 
Q for f2 = if one defines its amplitude reflection and transmission coefficients by 
cos g and sin g respectively while <fi denotes the phase difference between the reflected 
and transmitted fields . 

The parametric frequency conversion by a strong pump field (of frequency uo) in 

material can also be represented by an interaction Hamiltonian of the form Q 8 
with fi = 0. Here a and b are the annihilation operators for the signal (of frequency 
bj a ) and the idler (of frequency tot) respectively, g is a coupling constant that depends 
on the amplitude of the pump mode and (j) = Acot where Aw = lo + lo^ — uj a . We should 
note, however, that this Hamiltonian does not support parametric down conversion. 

Polarizing elements such as half- and quarter wave plates also can act as 
entangling devices. Quantum mechanically, polarized light is represented by a pair 
of orthogonal polarization modes (described by boson mode operators a, &), or as 
points on the Poincare sphere. The effect of a polarizing element on the field is a 
SU(2) transformation of the mode operators which corresponds to rotations on the 
Poincare sphere. The transformations are generated by Hamiltonians of the form £Q). 
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Finally, following the work of Wineland et al |U] , we consider a single laser cooled 
ion confined in a two dimensional harmonic trap. The internal and motional degrees 
of freedom of the ion can be coupled by applying two classical laser beams. If a and b 
represent the two oscillatory modes of the ion's quantized motion and <j> denotes the 
difference in phase between the two applied fields, then, under certain conditions 10 
the Hamiltonian for the ion's motion will be of the form Q in the interaction picture. 

The present work is also relevant in the context of parallel developments in the 
field of optical vortices. An optical vortex of order I centered at the origin (r = 0) 
has a field distribution of the form F(r) exp(il<j>). The distribution is such that the 
field intensity tends to zero as r — * whereas the phase shift in one cycle around the 
origin is 2ttI where I is an integer. The azimuthal mode index I has a physical meaning 
in that the vortex carries an orbital angular momentum of Yh per photon This 
angular momentum can be imparted to microscopic particles in order to manipulate 
them optically [121 1T3| . In recent years, this understanding has led to considerable 
interest in the generation and study of optical vortices both in free space ^31 an d in 
guided media [T51ITB] . 

A physically realizable field distribution that contains optical vortices is a higher- 
order Laguerre-Gaussian (LG) beam whose waist-plane field amplitude is given by 

E2 



u™(x,y,«>) = ,/XlJ)^ e -«(™-) (r V2/^)l m -"lLl" l -"l(2r 2 /^) e - 2 /- 2 (5) 

where r 2 = x 2 +y 2 , 9 = arctan(y/x), u is the beam waist, p = min(n, to) and L l p (x) is a 
generalized Laguerre polynomial. LG beams can be produced directly from a laser ^2] . 
In fact, in a hydrodynamic formulation of laser beam dynamics in terms of LG modes, 
vortices were found to occur in transverse laser patterns ^5] GO]- Usually however, 
LG beams are produced by the conversion or combination of Hermite- Gaussian {HG) 
beams that are emitted by most laser cavities. This is made possible because of the 
fact that any LG mode can be expressed in terms of HG modes. The waist-plane 
amplitude of the HG modes has the form 

u"m(x,V,u) = $ n (x,Uj)$ m (y,U)) (6) 

where 

$ n {x,w)=( ^ , 1 H n (V2x/w)exp(-x 2 /w 2 ) (7) 

and H n (x) is a Hermite polynomial. The decomposition of a LG mode in terms of 
HG modes is given as |T3] 

m+n 

Un%( x iV> u; )= ^i k b(n,m,k)u'^ ( l n ^ k:k {x 1 y,u) (8a) 

k=0 



(n + m)\k\ 1 d v „, , , , 

Kn, m ,k) =^i_ i ^__ jf[ (i-t ) »(i + tr]j« m 

The vortices as discussed above appear on the transverse amplitude profile of 
classical wave fields. Vortices can also occur in the configuration space representation 
of quantum systems of matter or radiation. Since the HG modes are also the energy 
eigenfunctions of a quantum oscillator, quantum vortices should arise in the study 
of wave packets of a quantum system that could be a two-dimensional harmonic 



Generation and evolution of quantum vortex states 



4 



oscillator like an ion in a two-dimensional trap. For a two-mode radiation field 
characterized by the annihilation operators a, b, and represented by a state vector 
\tp), the quantum vortex will appear in the quadrature distribution |(x, ylV')! 2 where 
\x,y) is the eigenvector of (a + a^)/\/2 and (b + &t)/\/2. Quadrature distributions 
can be measured by a homodyne method 21 . Vortices of matter will appear in the 
configuration space probability distribution. Recently it has been shown that the HG 
and LG modes are unitarily related[22] and the Poincare spherej2II] representing LG 
beams has an underlying SU(2) structure [21] • The Hamiltonian (f ) is therefore ideally 
suited to explore the possibility of generating quantum vortices. 

The objective and the plan of the paper are as follows. In section[2]we obtain the 
state vector and the wave function of a two-mode system which is initially in a Fock 
state and is acted upon by the Hamiltonian JJJ. We show how the two modes get 
entangled by evolution and under certain conditions evolve into a vortex state. The 
degree of entanglement between the modes is measured by finding the dependence of 
the von Neumann entropy on the system parameters. In section [21 the above analysis 
is carried out when the two-mode system is initially in a state that can be obtained 
from a Fock state via a unitary transformation. As an example, a quantum vortex 
is used as the initial state. We also find conditions for the revival and the charge 
conjugation of the vortex. In section 0] we consider the structure of the reduced state 
of each mode. The paper ends with concluding remarks in section 



2. Generation of quantum entanglement and creation of a quantum 
vortex using an initial two-mode Fock state 

2.1. Evolution of the state vector 

Let us consider the evolution of a two-mode Fock state \N—j, j) when the Hamiltonian 
is given by Q and the total number (N) of photons in the two modes is constant. The 
resulting state \ipNj (t)) = U (t)\N—j, j) can be obtained by the use of the disentangling 
theorem. In what follows, we use a different method. We write \N — j,j) as 

|W _ 3 ,, ) = fi^i |0 ,0, (9, 
and define a new pair of operators 

where U(t) = exp(-iHt) is the time evolution operator. Then \ipNj(t)) can be written 
in a compact form as 

lfcj(t)) = MM p) 

VN - f.f. 

Note that the state at time t is obtained by using the operators evaluated at time —t. 
The explicit expressions for a(t) and b(t) can be obtained by solving the Heisenberg 
equations for the operators. We get 

*D =v ® =v G) 
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where V = {i>«} is a 2 x 2 unitary matrix. Setting a — ^Q 2 + g 2 and £1 = crcosG, 
the matrix elements are written as 

vn = cos at — i cos G sin erf, V12 = — ie 1 ^ sin © sin erf, 

V21 = — ie _l ^smGsmcrf, V22 = cos at + i cos © sin at. (13) 

Note that 

Vll=V22, Ul2 = -«21 and l w 2l| 2 + 1^22 1 2 = 1- (14) 

Substitution in (fTT|> followed by binomial expansion and the use of © yields 

N-j j 

|V»JVj(*)) = X! E 6 ™™'^ - (m + n),m + n) (15) 

m— n.— 

where 

b mn=( N ~ j )( j )( N N Y /2 ( + ) 1/2 (^l) JV - J '- m (« 21 )" l («12)^"(« 22 )".(16) 

\ m J \n / \J\ — j J \m + nj 

The two modes in the state \ipNj) are entangled in the sense that the above double sum 
cannot be reduced to the product of two single-mode summations. 

It is instructive to briefly mention the case when the two modes are initially in a 
Glauber coherent state \a,/3). Since the Hamiltonian conserves photon numbers, 
the state at time t will also be a coherent state: 

U(t)\a,f3) = \a(t),f3(t)). (17) 

Applying Ijl2|l on \a,f3), we immediately obtain 

^-Vff). (18) 



Furthermore, the unitarity of V ensures that 

\a(t)\ 2 + \f3(t)\ 2 = H 2 + \(3\ 2 . (19) 

Thus no entanglement occurs if each input mode is in a coherent state. 

In what follows, we exploit coherent states as generating functions of number 
states to reduce the double sum in (fTK|) to a single sum. Expanding both sides of (|T7jl 
in number states and recalling that |a(i)| 2 + |/?(i)| 2 = |a| 2 + |/3| 2 , we get the relation 

fyin on n-P+'J 

m n v p q v 



where r = f3/a and 

Wr) = («n + v 1 2tT(v 21 + V22TY = E ^di k H Pq {r)\ T ^. (21) 

k=0 

Substituting in (|20|l and equating the coefficient of a N ~^ ^ , one gets 

N 

\i, N] (t)) = U(t)\N-j,j) = £tf$|JV - q,q) (22) 

9=0 

where 

1 



^(9) _ 
_ j! 



Some useful properties of |C^| 2 are derived in appendix A. 
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2. 2. The wave junction- a coherent superposition of vortex states 

The corresponding wave function in configuration space is obtained as follows. Using 
the relation 

P -y 2 / 2 H (v) 2 , s 2 

(y\q) = A H q (y) = (-l)'e* d^e^ , (24) 

and the corresponding expression for (x\N — q), we obtain 

tp Nj (x,y,t) = (x,y\U(t)\N - j,j) 

_ e-^+y 2 )/^ Mq) H N ^ q {x)H q {y) 

The wave function has a more appealing form in polar coordinates as shown below. 
Writing x — r cos 6, y = r sin 8, and defining 

7±( T ) = vn + vi 2 t ± i(v 2 i + v 2 2t), (26) 
we get (see appendix B) 

N 

ip Nj (x,y,t) = y^by]ttjy_ W)rt (r,g) (27) 

where 



{n) _ 1 / (iV-i)!i! (j) , , 



Cnu(t) = 7+ (r) w - 7 _(r)" 

«mn(r,9) =u^(x,y,V2). (28) 

Recall that for m ^= n, u mn (r, 9) represents a vortex of order \m— n\ and charge m — n 
embedded in a Gaussian host beam of waist ui = y/2. Thus for odd values of N , the 
wave function ipNj(x,y,t) becomes a coherent superposition of vortex states, whereas 
for even values of N , the superposition will also contain a state ( corresponding to 
n = N/2) that does not have a vortex character. |25j. 



2. 3. Creation of a single quantum vortex 

In this section we will derive conditions for the creation of a single quantum vortex. 
We reiterate that for light fields, the vortex will appear in the quadrature distribution 
whereas for other systems it will be in the probability distribution in configuration 
space. 

The initial two-mode Fock state can evolve into a single vortex state when the 
summation in (|27|l collapses into a single term. This happens whenever 7 + (0) or 7_ (0) 
is zero. It is easy to show that |t>2i| 2 = 1/2 for both these cases. 
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If 7+(0) = 0, then vn + iv 2 \ — and taking the complex conjugate of this 
equation, v 22 + = 0. Then 7+(t) = 2w 12 r and 7_(t) = —2iv 2 i so that 
C^1(0) = {2v 12 ) N - n {-2iv 2 i) n j\5 N - , n and, finally 

^■(x,J/,f)J 7+( o)=o = 2 JV / 2 l ^ A ' U21 -^ J 12 u J - A r_,(r,0). (29) 

The condition 7+(0) = implies that Q, = — gsinc/> and cr cos erf = — g cos sin erf. We 
give two examples for which these conditions are satisfied. 

(i) Setting £1 = 0, = it and at = 7r/4, we get 

ipNj(x,y,t) = i 3 u j:N ^j(r,6). (30) 

From H22f> . one obtains the corresponding state vector 

N 

U Q \N- J , 3 )=J2D ( £]\N~q 7 q) (31) 
g=0 



where 



and 



U = exp[^(a t 6 + a6 t )] (32) 



n(9) _ 1 

Nj ~ ( - / AT : /Jn=O,0=7r,(Tt=7r/4 

- ■ / mJ^A WO- + <r) w -«(l - iry) T ^ . (33) 



V 2 N (N-q)\q\j\ 
(ii) Setting fi = g, <j> = —tt/2 and erf = 7r/2, we get 

ip Nj (x,y,t) = (-iy +N u ]tN ^(r,0). (34) 
The operator form of the corresponding state vector is given by |26| 

exp[^={i(aib-atf) - (a 1 a - tfb)}]\N - (35) 
2v 2 

Following a similar analysis for 7_ (0) = 0, one obtains 

ipNj(x,y,t)\ 7 _( )=o = 2 N/2 i N ~ 3 v 21 ~ 3 <v{2UN-j,j{r,6). (36) 

The condition 7-(0) = yields £1 = gsm<f> and a cos at — g cos c6 sin erf. These 
two conditions are satisfied , for example, when f2 = cf> = and erf = 7r/4. The 
corresponding wave function is the complex conjugate of J30J. 

We end this section by noting that the above conditions can be physically realized 
for a given entangling device. We give an example in the context of a frequency 
converter. Suppose the signal (of frequency w„) and the idler (of frequency LUb) are 
initially in Fock states and the converter is pumped at the difference frequency uj a — tOb- 
Replacing f by L /c, where L is the length of the non- linear medium and c is the speed 
of light, one can adjust the pump amplitude such that gL/c — 7r/4. This setup 
corresponds to il = <j) = and gt = 7r/4. In this case the quadrature distribution of 
the output state will be a single quantum vortex as mentioned above. 
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2.4- Entanglement of the two modes 

Initially the two modes are not entangled as the state vector \N — is the direct 
product of the state vectors for each mode. In configuration space, this would imply 
that ipNj(x,y,0) is separable in x and y as indeed it is. Furthermore, as the time 
dependence arises solely in iiy which vary as cos at or sin at, the initial state is revived 
whenever at = kir where k is an integer. For even values of fc, the revival is exact 
whereas for odd values of k, it is within an overall factor of (— ) N . At other times, the 
two modes are entangled as is evident in the expression C (|15|l or (|22|l 'l for the state 
vector and the expression f i|25|l or $27$ ) for the corresponding wave function. 

2.5. Degree of entanglement 

Note that the two-mode system \ip^j(t)} is in a pure state whereas the reduced state 
of each mode, determined by a partial trace operation, will be a mixed state. The 
reduced density operators of modes 'a' and 'b' are given respectively by 

JV 

= Tr b \i> Nj )(i; Nj \ = J2 \C$\ 2 \q)(q\ (37a) 

9=0 
N 

P% = Tr a \i> Nj ){^ Nj \ = ]T \C { N N r^\\)( q \ (376) 

9=0 

The corresponding von Neumann entropies and provide a measure of the 
degree of entanglement between the two modes: 

N 
9=0 

= -El^ ?) I 2 log|C^| 2 (386) 

9=0 

By virtue of relations (|A.4(1 , we get 

£>Nj\ \v2i\ 2 ^l-B = £>Nj\ \v2i\ 2 ^R = ^N,N-j\ \v2i\ 2 ^R- (39) 

Changing the summation index from q to N — q in the expression for , one obtains 

^n] = ^Nj- Thus the symmetry relations hold good for S 1 ^ as well. These 
observations hold for any bipartite system in a pure state. 

It is remarkable that for a given value of N, j and q, the dynamics of |C^ | 2 
depends on |i>2i| 2 = sin 2 0sin 2 at only (see appendix A). This important observation 
implies that (a) the entropy S^j and the reduced density operator are independent 
of <p and (b) are symmetric with respect to the interchange of <d and at. In Figure 1, 
we plot S^j as a function of |t>2i| 2 for N = 4 and j — 0, 1, 2. 

Trivially, for |i;2i| 2 = 0, the initial pure state \N — j,j) either does not evolve or 
is fully revived and the entropy of the reduced state is zero. For |w2i| 2 = 1, the initial 
state swaps the photon numbers in the two modes and becomes \ j, N—j) which is also a 
pure state. For all other values of |^2i| 2 , the initially pure state becomes a mixed state 
and the entropy of the reduced state becomes non-zero. Recall that for |i>2i| 2 = 1/2 
and N — j ^ j, the quantum state becomes a vortex. Thus a quantum vortex is indeed 
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Figure 1: Plot of S$ as a function of \v 2 i | 2 for N = 4 and j = 0, 1, 2. 
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Figure 2: Plot of as a function of j for \v2i\ 2 = 1/2 and a given total number of 
photons N. 



an entangled state. To quantify the degree of entanglement for a vortex state, we plot 
Sffj as a function of j for \v2i | 2 = 1/2 and a given total number of photons N (see Fig. 
2). It is clear that the entropy of the state without a vortex (j = N/2) is less than 
the entropy of the neighboring (j ~ N/2) vortex states (N — j ^ j). This reduction 
in entropy can be attributed to the symmetry of the j — N/2 state and traced to 
the highly oscillatory nature of the Jacobi polynomial appearing in Eq. i|A.5|) . For a 
given value of N, the minimum in the entropy of a vortex state occurs for j = 0, N 
in which case |C^| 2 is a binomial distribution (see appendix A). Interestingly, for 
j = 0, N, the vortex state will have the maximum allowed order (N). Thus the vortex 
state of maximum order will have minimum entropy which is counter-intuitive. One 
would have expected that the more twists the phase of the state has, more energetic 
and more entropic it would be. Note that the symmetry of about |«2i| 2 = 1/2 in 
Fig. 1 and about j — N/2 in Fig. 2 is contained in the relations (|39|l . Note also that 
the vorticity or non-vorticity of the state of lowest entropy depends on the value of N 
(see Fig. 3). 

We end this section by comparing the entropy values in Figs. 1-3 with log 2 (A r +l), 
the maximum entropy possible for a given N with an entirely mixed state. For N = 4, 
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Figure 3: Plot of S^j as a function of N for \v 21 \ 2 = 1/2 and 
j = 0,1, N/2. 



10 and 100, \og 2 (N + 1) has the values 2.32193, 3.45943 and 6.65821 respectively. 

3. Evolution of an initial vortex state 

3.1. Evolution of the state vector 

Let us assume that the two modes are initially in a quantum vortex state as in (|31|) . 
Then the state vector at time t will be given by 

$ Nj (t))=U(t)U \N-j,j) (40) 

Proceeding as in section |2~T1 we define 

(4i) 

and obtain 

/a(t)\ = fa(0)\ 

\kt)) Ui 

Note, however, that in this case 



u (t) r v U(o)J- (42) 



Thus 

Q "*(»)' ™ ,44) 

The action and the effect of the unitary operator Uq are now clear. Uq transforms 
the two-mode Fock state into a different initial state before its time evolution begins 
and thus Uq can be regarded as the operator for initial state preparation. The effect 
of Uq is contained in the unitary matrix W. As a result, the overall unitary evolution 
matrix changes from V to V = VW. 

In the present case, /To, as given by l|32|l . prepares a quantum vortex state as the 
initial state and the corresponding expression for W is given as in (|43() . A different 
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expression for Uq will generate an initial state that is different from a quantum vortex. 
Yet for all these initial states the dynamics is essentially solved once the corresponding 
dynamics for a two-mode Fock state is worked out as in section IzTTl In each case, one 
need only calculate the matrix W and replace V by V. In this sense, our theory not 
only provides a unified approach to entanglement through a generic Hamiltonian but 
also promises wide applicability to a variety of initial states. 

In the present case, the matrix elements of V = {"y} are obtained easily as 

"11 = (vu + ivu)/V2, V12 = ("12 + it>ii)/v2, 
"21 = ("21 + 2^22)/ V2, V 2 2 = ("22 + iV2i)/v2. (45) 

Using 114(1 . one can also show that 

"11 = "227 "12 = -"21 and l"2i| 2 + l"22| 2 = 1- (46) 
It is now trivial to obtain the wave vector and the wave function by borrowing 
the corresponding results from the previous section. We simply replace by for 
i,j = 1, 2 and for the sake of clarity and comparison, use the same nomenclature for 
the new expressions except for a~ (tilde) over them. Thus 

JV 

i«<)> = E^>-^> ( 47 ) 

9=0 



where 



1 



r K — 



_(N- q )\ q \\ 9 ^-™(^o (48a) 

W T ) = («n + "i2t) p (w 2 i + "22T) 9 (486) 
Furthermore, the results of appendix A can be used to write 

1 2 \ 1-3 



\C ( ^ = (N - j)l(N - q )l q \Ur\l - \i2i\T [j^^) l/g(|S2i| 2 )| 2 (49) 

(50) 



8q,j, |"2l| 2 ^0, 
$q,N~j, |"2l| 2 ^l- 



and 



\ C n]\ 2 \\v2i\ 2 ~*1-R - lAvj q) \ 2 \\v 2 i\ 2 ^R - \ C n]n -]\ 2 \\v2i\ 2 ^B.- ( 51 ) 



3.2. The wave function 

The corresponding wave function in configuration space can be read off from Eq l(27|l . 
We get 

N 

i>Nj{x,y,t) = y^y$]uN- n ,n{r,0) (52) 

where 




(N-jy.ji z U ) , , , , 

(7V-n)!n!2^ Cjv » (U) {56a) 

UW=7 + W"-"7-(t)" (536) 
7 ± (t) = in + vi 2 T ± i(v2i + V22t). (53c) 
Thus a quantum vortex state evolves into a superposition of vortex states under 
the action of the Hamiltonian ||IJ. 
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3.3. Revival and charge conjugation 

It can be shown that if 7+(0) or 7_(0) is zero, then the summation in (|52|l reduces to 
a single term. Specifically, if 7+(0) = 0, then Imi^i = Im«22 = and 

4>Nj(x,y,t) = {iv) : >v* N ~ : 'u : j. N -j{r,e) (54) 

with v = Ret>22 + iReuai- The above conditions are satisfied for the following cases: 
(a) sin at = for arbitrary values of O and 0. This includes the initial state (t = 0) 
and the state upon revival (at = ir). (b) sin0 = sin0 = 1 for arbitrary time. In this 
case the initial vortex state becomes an eigenstate of the corresponding Hamiltonian. 
On the other hand, if 7-(0) = 0, then Re«2i = Re«22 = and 

$ Nj (x,y,t) = (v¥(-iv*) N -j UN _ jd (r,6) (55) 

with v = Imw22 +ilmu2i. Note that uat_ (r, 0) — u* N ^(r,9) and thus 7-(0) = is 
the condition for 'charge conjugation' or 'helicity reversal' of the initial vortex state. 
This condition is fulfilled whenever sin O sin = and cosci = 0. 



3.4- Degree of Entanglement in the superposition state \bty) 

The reduced density operator of mode 'a' and the corresponding von Neumann entropy 
are given respectively by 

N 

p$. = Tr b $ Nj )$Ni\ = £ Pn]\ 2 \q)(q\ (56a) 

9=0 

N 

41= -EieSg| a log|C<g| a (566) 

q=0 

It is clear that S^j will depend on |w2i| 2 in exactly the same way as S^j does 
on |^2i | 2 except that the form of |v2i| 2 as a function of 9, <fi and at is quite different 
from | V21 1 2 - Notably, |u2i| 2 depends on <j> while |t>2i| 2 does not. Explicitly, 

|i>2i| 2 = — — sin O sin at(cos<f> cos at — sin <f) cos sin at) . (57) 
3. 5. A case of constant entropy 

Note that if 9 = cf> = tt/2 or 6 = 0, then |u 2 i| 2 = 1/2 so that |C^]| 2 and the entropy 

Sftj will not evolve with time. The underlying reason is as follows. 
Using the expressions @ for the SU(2) generators, we obtain, 

J 3 \N- q,q) = ?t^l\N-q,q) (58) 

Noting that Uo, as given bv I32fl . can be written as Uq — exp(i7rJi/2) and using the 
relation J 2 — UqJzUq, we also get 

J2U \N-j,j) = ^^-U \N-j,j) (59) 

For = <f) = tt/2, the Hamiltonian reduces to H = —2gJ^ for which Uo\N — 
becomes an eigenstate by virtue of f5§|l . 

The condition = corresponds to g — and the Hamiltonian reduces to 20,,]%. 
From lEU, EH) and gZJ, one then immediately obtains 6%) = exp(-ffit[iV-2g])Z)^ 
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Figure 4: Time evolution of an initial vortex state for O = even though the entropy 



Sj^- remains constant. Shown here are the contour plots of the absolute square (top row) 



J Nj 

and the phase (bottom row) of ipNj(x,y,t) as functions of x and y at different times with 
bffj as given in 160H . Here, N = 4, j = and cos 2 fit has the values (a) 1.0, (b) 0.9, (c) 
0.5 and (d) 0.0. The horizontal and the vertical axes refer to the x and y coordinates 
respectively. For the phase plots, we have used the convention that the phase ranges from 
— it to 7r. Note that for cos 2 fit = the state becomes the complex conjugate of the initial 
state as the direction of phase change is reversed. 



so that |C^]| 2 , and consequently, the entropy S^j become independent of time. The 
corresponding wave function is given by (|52|l where the coefficients have the value 



l(n) 



-1/2 



N 



-1/2 



H ) iv (sin fit) N - n+j (cos flty 1 -! f$! (cos 2 fit) . (60) 



(n), 



It is interesting that although IC^I 2 and the entropy remain constant for = 0, 
the initial vortex state will continue to evolve with time as shown in Figure 4 [?7 . For 
TV = 4 and j = 0, the initial state at t = (Figure 4(a)) corresponding to cos 2 fit = 1.0 
is a vortex of order 4 and charge -4 as given by equations (|30|l and (|31Jl . Recall that 
fl = (7 cos 0. Thus = corresponds to fl = a. Furthermore, if cos 2 fit = 0.0, then, 
with = 0, the conditions for charge conjugation as given below equation l|55|) are 
satisfied and we obtain the complex conjugate of the initial vortex. For cos 2 fit = 0.5, 
it is more convenient to use cartesian co-ordinates. Using the expression for as 
given above and the configuration space representation of number states as given by 
J23J, one can use the summation theorem for Hermite polynomials 29 to obtain 



^Ao{x 1 y,t)\ nt= 



-(x 2 +y 2 )/2 



Vtt2 4 4! 

e -(x 2 +y 2 )/2 



-H 4 



y 



V2 



(x-y) 4 + 6(x-y) 2 -3] 



(61a) 



(616) 



Thus, the wave function is a Gaussian modulated by a Hermite polynomial. Clearly, 
its value is real and, therefore, its phase is either zero or tt depending respectively 
on whether the wave function is > or negative. It is easy to show that the wave 
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function vanishes whenever (x — y) 2 — 3 ± \f&. Finally, it may be of some interest to 
realize that the wave function corresponds to a SU(2) coherent state —\t,N) in the 
Schwinger representation with r = — 1 and N = 4. 

4. Structure of the reduced state 

In order to determine the structure of the reduced state for each mode, we first consider 
the correlation function in the rr-space of mode 'a' given by (x\p^\y). A classical 
analog of this function is the mutual coherence function of a partially coherent source 
28 . Thus the process of reduction of a pure two-mode state into a mixed state by a 
partial trace operation over one mode amounts to loss of coherence and information. 
We can also define the spatial coherence function 7^j(0 for the reduced state by 

7$(0 = f(x\ P { ±x + l)dx (62) 



When the system is initially in a two-mode Fock state, one obtains 

N 1^(9)12 

top®\v) = E y^=e-^ + y 2 ^H q (x)H q (y). (63) 

The corresponding expression for 7^j(i) is obtained by evaluating the standard 
integral [H| in (52). We get 



N 



" (a) (0=El C Sl 2e " i2/ %^/2) (64) 



q=Q 



where L q (x) = L q (x) is a Laguerre polynomial. Note that only one term survives in 
the summations over q by virtue of (|A.3|> whenever |t>2i| 2 = or 1. Thus 

e -(*W)/2 f H ^ {V \ hi| 2 -0, 
(x\M\y) = - o ' _ _ (65) 



and 



-i(N-j)\ 



V21V -» 1- 



1 P(l) = e - l2 /U - (66) 
I L N ^(P/2), \v 21 \ 2 ^l. 

Furthermore, one can use Eq. IjATf) and the definition (|37af) to get 

{x\pN]\y)\\v 2 i\ 2 ^l-R = ( x \PN,N-j\y)\\v 2 i\*^R ( 67a ) 

7^j(0J|t>2i| 2 -i--R =7^,jv-jJ|t>2i| 2 ^R ( 67& ) 
When the system is initially in a vortex state, Eqs. H63I676(I are still valid provided 
that |w2i| 2 is replaced by |"D2i| 2 - In Fig. 5 we present contour plots of the correlation 
function as a function of x and y for a set of values of |«2i| 2 when N = 8 and j = 
(top row), j = 4 (bottom row). The intricate patterns for |«2i| 2 = and 1 can be 
explained by using Eq. (|65|l . Furthermore, the identical nature of patterns for N = 8, 
j = 4 on cither side of |«2i| 2 = 1/2 can be attributed to the property (|67afl . Finally 
in Fig. 6 we plot 7^(0 as a function of I and |«2i| 2 when N = 4 and j = 0, 2. The 
patterns for |i>2i| 2 = and 1 follow from Eq. I|66|) and the symmetry of the plot for 
j = 2 about \v 21 \ 2 = 1/2 follows from (j676j) . 
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Figure 5: Contour plots of the correlation function as a function of x and y for different 
values of |«2i| 2 - The parameters are as follows: N = 8; j = (top row), j = 4 (bottom 
row); |v2i| 2 has the values (a) 1.0, (b) 0.9, (c) 0.5, (d) 0.1 and (e) 0.0. The horizontal and 
the vertical axes refer to the x and y coordinates respectively. 
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Figure 6: Contour plot of 7^"j(0 as a function of / and |i>2i| 2 when N = 4 and j = 
(left), j = 2 (right). 

5. Conclusion 

In conclusion, we have studied, in a general way, entanglement produced in a two- 
mode bosonic system by linear SU(2) transformations leading to the generation and 
evolution of quantum vortex states. The linear SU(2) transformations are generated by 
evolving the system under the action of a generic Hamiltonian that mimics a variety of 
entanglement mechanisms. We have demonstrated that these transformations produce 
a coherent superposition of quantum vortices in general, and a single quantum vortex 
under certain conditions. Furthermore, as one would expect, a vortex state is found 
to be an entagled state. When the system is a light field, the vortex will appear in the 
quadrature distribution that can be measured by a homodyne method .21 . Explicit 
analytical results were obtained when the system was initially either in a Fock state 
or in a quantum vortex state. In the latter case, we have also found conditions for 
its revival and charge conjugation. A simple recipe was provided to accommodate all 
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other cases for which the initial state can be reached from a Fock state by a unitary 
transformation. Thus we not only provide a unified approach to entanglement through 
a generic Hamiltonian but also predict wide applicability of our results to a variety of 
initial states. 

The ideas developed in this paper can be applied not only to light fields but 
also to matter waves such as the Bose Einstein condensates (BEC). In recent years, 
the BEC has proved to be an excellent laboratory for studying (both bipartite and 
many-particle) entanglement j-'iOl 1311 1321 1331 IHl] . The entanglement of the modes as 
well as the entanglement of the atoms in a BEC have been considered. We mention 
parenthetically that our work is relevant in the former case. It is also well known 
that several mechanisms exist for the generation of vortices in a BEC 35 , 3BJ |23 EH| • 
Additionally, Whyte et al. jSH] have used the similarity between BECs and laser 
light to propose a method for generating Hermite-Gaussian type modes in a so- 
called light pulse resonator. Thus it should indeed be possible to generate vortices 
in a two-component BEC by entangling the two modes of the BEC by linear SU(2) 
transformations via an entangling device such as a beam splitter. 
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Appendix A. Some useful properties of |Cj^| 2 

Using Leibniz' rule for the j-th derivative of a product and the relations lfH|) . one 
obtains 

|cg| 2 = (N-j)\(N - ?)!g!C?!) -1 (l - \v 2 i\ 2 ) N (t^^p ) ? ' l/$(l«ii| a )| a (A- 1 ) 
with 



It is easy to show that 



|c7g|2= i r , , 2 i (A - 3) 

Next we derive some important symmetry properties of |Cjs^j| 2 . First we show that 

where < R < 1. The first relation is obtained from l)A.2(l by changing the summation 
index from ktoj — k and the second relation is proved by exploiting the non-negativity 
of the factorials in (|A.2() and changing the summation range accordingly. Using these 
two relations we immediately get 

\ C Ni?\\vii\ 3 ->l-R = \ C Nj 9 \ 2 \\v 21 \ 2 ^R = \ C N,N-j\ 2 \\v2i\ 2 ^R- ( A - 4 ) 



5a,i, |W21| 2 ^0, 
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Note that if \v 21 \ 2 = 1/2, then |u n | 2 = \v 22 \ 2 = \v 12 \ 2 = 1/2 as well. Additionally, 
if j = or AT, then \C^]\ 2 = 2- N ( N ) is a binomial distribution whereas if j = N/2, 
then 

\C$ 2 , N/2 \ 2 = (TV!)^[W2)!P^ 2 /2 -^- JV/2) (0)] 2 Q (A.5) 
where P n a ' l3 \x) is a Jacobi polynomial. 
Appendix B. Derivation of equation (|27[) 

Substituting the expression I12,'S[1 for C^j in (|25|l and performing the summation over 
q before differentiation with respect to r, we get 

M*,V,t) = ^^|^e^^ 2 )/ 2 [^)i^( T ) e -(^+^)]_ (B.l) 
where 

i(r) = (u u + ui 2 r)<9 x + (u 2 i + t^r)^. (B.2) 

We introduce z — x + iy so that x 2 + y 2 = zz* and A(t) = 7 + (r)9 z + 7_(r)9 z . with 
7± given by (|26|) . Next we expand A(r) iV binomially and then use the relation 

f {-l) n m\e- zz " z n - rn L™- m {zz*), m < n, 
dl n d 7 z l ,e- zz * = { " (B.3) 

1 (-l) m n\e- zz * z* m -' n L™- n (zz*), n < m. 

to evaluate A(t) n e~( x2+y2 \ Collecting all the terms we finally obtain l|27() . 
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